Manuscript for Orbital2001 



Magnetism and Superconductivity in a Two-band Hubbard Model in 
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We study a two-band Hubbard model using the dynamical mean-field theory combined with 
the exact diagonalization method. At the electron density n = 2, a transition from a band- 
insulator to a correlated semimetal occurs when the on-site Coulomb interaction U is varied for 
a fixed value of the charge-transfer energy A. At low temperature, the correlated semimetal 
shows ferromagnetism or superconductivity. With increasing doping \n — 2|, the ferromagnetic 
transition temperature rapidly decreases and finally becomes zero at a critical value of n. The 
second-order phase transition occurs at high temperature, while a phase separation of ferromag- 
netic and paramagnetic states takes place at low temperature. The superconducting transition 
temperature gradually decreases and finally becomes zero near n = 1 (n = 3) where the system 
is Mott insulator which shows antiferromagnetism at low temperature. 

KEYWORDS: two-band Hubbard model, metal-insulator transition , ferromagnetism, antiferromagnetism, super- 
conductivity, dynamical mean-field theory, infinite dimensions, exact diagonalization method 
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§1. Introduction 

The two-band Hubbard model has been extensively 
studied as a simple model simulating transition-metal 
compounds and high-T c superconductors, which is char- 
acterized by two parameters: the on-site Coulomb inter- 
action U at the ei-orbital and the charge-transfer energy 
A between the d- and p-orbitals. Recently, several au- 
thors have studied the model using the dynamical mean- 
field theory (DMFT) which becomes exact in infinite spa- 
tial dimensions. 1_7 ^ The Mott metal- insulator transition 
was found to occur at n = 1 (or n = 3), 1:3 ' 5 ~ 7 ) where n is 
the total electron number per unit cell and given by the 
sum of p- and <i-electron numbers: n = n p + n<j. Anti- 
ferromagnetism 4 ) and superconductivity 1,2 ' 5 -* were also 
observed in this model. However, the ferromagnetism 
was not discussed there. 

Recent findings of the high-temperature ferromag- 
netism in doped CaBg and the high-temperature super- 
conductivity in MgB2 have simulated a renewed inter- 
est in strong electron correlation in the multi-band sys- 
tems. Therefore, we study the ferromagnetism and the 
superconductivity in the two-band Hubbard model us- 
ing the DMFT combined with the exact diagonalization 
(ED) method. 2, 5 ^ We especially focus on the electronic 
states near n = 2 where the system is semimetal or band- 
insulator in the non-interacting case. 

§2. Formulation 

The two-band Hubbard model is given by 
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where d\ a and Pj a are creation operators for an electron 
with spin a in the d- orbital at site i and in the p- orbital 
at site j, respectively. The charge-transfer energy A is 
defined by A = e p — e<j. In eq.(2.1), we assume that p- 
and d-orbitals are on different sub-lattices of a bipartite 
lattice. 

In the limit of infinite dimensions, the self-energy 
becomes purely site-diagonal and the DMFT becomes 
exact. The local Green's function for the d-electron, 
D a (r — t') = —{Td a {T)d\{T r )) , can be given by the im- 
purity Green's function of an effective single impurity 
Anderson model, 

H And = e f „fiu + uf\hf\h 

+ E e kJ ka c k(J + £ V ka (p a c ka + 4/ CT )(2.2) 

k , a k,a 

where Ef a is the impurity level and e ka are ener- 
gies of conduction electrons hybridized with the im- 
purity by V ka . In the model eq. (2.2), the non- 
interacting impurity Green's function, C?o CT (iw n ) = 



, includes effects of the in- 



k iuj n — eka 

teraction at all the sites except the impurity site and 
is determined self-consistently so as to satisfy the self- 
consistency equation. 

For simplicity, the calculations in this paper are re- 
stricted to the Bethe lattice with the connectivity z and 
the hopping tij = and we set t p d = 1. In the 
limit z — co, the self-consistency equations for the lo- 
cal Green's functions are given by 1 ) 



Goa(Su n ) 1 
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= iu n + fj,- e d - t^ d P a (kL> n ), 
= iw„ + pi - e p -t* d D a (iu} n ), 



(2.3) 



2 



Yoshiaki Ono and Kazuhiro Sano 



u 10 



Insulator | 




n=2 T=0 













,---'8' 






A=6 T=0 



10 



Fig. 1. The phase boundary separating the metallic and insulating 
regimes as a function of U and A at n = 2 and T = 0. The inset 
shows the chemical potential as functions of n for U = 7,8, 9, 10 
at A = 6 and T = calculated from the DMFT (ED method) 
with N s = 8. 

where P a (iui) is the local Green's function for the p- 
electron and \i is the chemical potential. 

To solve the impurity Anderson model eq.(2.2) we use 
the exact diagonalization of a finite-size cluster with Eko 
and V ka for k = 2,3,. ..,N S (ED method). 2) First, we 
assume 4N S — 2 parameters: {sf a , £ka, Vko-}oid- Next, 
we diagonalize the finite-size cluster with these param- 
eters, and calculate the local (impurity) Green's func- 
tion Da-(ioj n ). Finally, we calculate Goo-fan) through 
the self-consistency equations (2.3), and determine the 
new parameters {sf a , Ska, Vk a }ncw so as to make the non- 
interacting Green's function with finite N s , £,^f(ia; n ), as 
close to Go<r(iu n ) as possible, namely, so as to minimize 
X 2 defined by x 2 = E w „ \Goa(iu n )~ 
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where Q^i^n) = (iw„ - e fa - ^ k=2 lu . 
process is iterated until the solutions converge 
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§3. Metal-Insulator Transition at n = 2 

First, we consider the paramagnetic state at zero tem- 
perature for n = rid + np = 2. In the non-interacting case 
U = 0, the system is a band-insulator for A ^ while it 
is a semimetal for A = 0. Within the restricted Hartree- 
Fock approximation (HFA), the energy gap is given by 
A - Then the system is metallic for U = 2 A, oth- 
erwise it is insulating. We note that rid = n p = 1 for 
U = 2A due to the particle-hole symmetry. 

In the inset in Fig. 1, we show the chemical potential 
H as functions of n for U = 7, 8, 9, 10 at A = 6 and 
T = calculated from the ED method with the system 
size N s = 8. In this calculation, the solution is restricted 
to the paramagntic state with £f a — Sf, eta = £k and 
Vka = Vk for k = 2,3, ...,N S . When U increases from 
U = 0, the discontinuity in the chemical potential at 
n = 2 decreases and finally becomes zero at a critical 
value, where a transition from the band-insulator to the 
correlated semimetal occurs. 1 ) The critical values for 
the metal- insulator transition are plotted in Fig. 1. In 
contrast to the restricted HFA, the metallic state is found 
in the wide parameter region due to a correlation effect 
considered in the DMFT. 
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Fig. 2. The magnetization for the d-electron M^, that for the p- 



electron M v and the total magnetization M 



functions of the temperature T at n = 2, obtained from the ED 
method (closed circles) for U = 8 and A = 4 with N s = 6 and 
from the Hartree-Fock approximation (dashed lines) for U = 2 
and A = 1. 
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Fig. 3. The total magnetization M as a function of the electron 
number n, obtained from the DMFT for U = 8 and A = 4 at 
T = 0.01 (closed circles), 0.025 (open circles), and from the HF 
approximation for U = 2 and A = 1 at T = 0.025 (dotted line), 
0.1 (dashed line). 

§4. Ferromagnetism 

At low temperature, the correlated semimetal men- 
tioned above becomes unstable compared to a ferromag- 
netic state. In Fig. 2 we plotted the magnetization for 
the d-electron Md, that for the p-electron M p and the 
total magnetization M = Md + M p as functions of the 
temperature T at n — 2 for U = 8 and A = 4 calcu- 
lated from the ED method with N s = 6. As shown in 
Fig. 2, Md and M p have opposite sign to each other. In 
the low temperature limit, both of Md and M p become 
constant, while the sum of them M becomes zero. The 
feature of the ferromagnetism from the DMFT is similar 
to that from the HFA as shown in Fig. 2. However, the 
transition temperature T c from the HFA is much higher 
than that from the DMFT (see also Fig. 4). 

Fig. 3 shows the total magnetization M as a function 
of the electron number n. When n decreases, M contin- 
uously becomes zero at a critical value of n for high tcm- 
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Fig. 4. The transition temperature T c for the ferromagnetism as 
a function of the electron number n, obtained from the DMFT 
(closed and open circles) for U = 2A = 8, 4, 2, and from the HF 
approximation (solid and dotted lines) for U = 2A = 2. 

peratures (see for T = 0.025), while it discontinuously 
becomes zero for low temperatures (see for T = 0.01). 
The similar properties are also observed within the HFA 
as shown in Fig. 3. At low temperatures, however, the 
HFA also predicts a metastable state where M decreases 
with increasing n and continuously becomes zero at a 
critical n (see for T = 0.025). By calculating the ther- 
modynamic potential, we find that the phase separation 
of the ferromagnetic state and the paramagnetic state oc- 
curs at the low temperatures as shown in Fig. 4. The pos- 
sible metastable state and the phase separation within 
the ED method will be reported elsewhere. 

Fig. 4 shows the transition temperature for the fer- 
romagnetism T c as a function of n for several values of 
U(= 2A). T c monotonically decreases with decreasing 
n. The closed circles show the second-order phase transi- 
tion, while the open circles show the discontinuous tran- 
sition as seen in Fig. 3. Within the HFA, the second- 
order phase transition occurs at the high temperature 
(solid line) , while the phase separation occurs at the low 
temperature (area between the dotted lines). 

§5. Superconductivity 

Finally, we discuss the superconductivity. 1,2,5 ' The 
on-site paring susceptibility x P is given by 1 ' 

X V = ^ f &tYJ?*K (t)*! (7-^(0)^(0)) 
= Tj2l^ 1/2 {I A}- 1 • A • a-^.y, (5.1) 

v.y' 

where [A] u y = t* d \P(w)\\x\oc]vy\P{iv')\, and \ a ] v y = 
tpd\P(S u )\ 2 5 v y ■ Here, xi oc is the local paring susceptibil- 
ity on a d-orbital given by 

rP rP rP rP 
[XiocUv = T 2 / dn / dr 2 / dr 3 / dne 1 ^-^ 
Jo Jo Jo Jo 

xe -'(r 3 -r 4 ) (rdT(ri)di(T2)d t (T3)d t (r4) ^, 2) 

To calculate \\ oc within the exact diagonalization 
method, we use a spectral representation of r.h.s. in 



0.02 



0. 



1 1 1 1 1 

U=8 A=4 








1 



Fig. 5. Transition temperatures for the ferromagnetism T c , the 
singlet superconductivity T c and the antiferromagnctism Tjv , ob- 
tained from the DMFT for U = 8 and A = 4. 

eq.(5.2) by inserting a complete set of eigenstates 
When the largest eigenvalue of A approaches unity, the 
paring susceptibility diverges. It signals the transition 
into the superconducting state from the normal state. 

In Fig. 5, the transition temperature for the singlet 
superconductivity T s is plotted as a function of n for 
U = 8 and A = 4, obtained from the ED method with 
N s = 4. 5 ) We also plotted the transition temperature 
for the ferromagnetism T c calculated in §4 together with 
that for the antiferromagnctism Tjv calculated from the 
ED method with N s = 6. 

§6. Discussion 

We have obtained a phase diagram including the 
metal-insulator transition, the ferromagnetism, the an- 
tiferromagnctism and the superconductivity using the 
DMFT combined with the ED method, for the two-band 
Hubbard model with typical values of the parameters U 
and A. More systematic calculations for various values 
of the parameters are now under the way. 

Effects of the hopping integrals between the p-p and 
the d-d orbitals, which are not considered in the present 
study but are not negligible in actual compounds, make 
important contribution to the transition temperatures. 
Such effects and the detailed analysis of the mechanism 
for the superconductivity will be reported in a future 
publication. 
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